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Abstract 

We propose a dipole modulation model for the Cosmic Microwave Background Radiation (CMBR) po¬ 
larization field. We show that the model leads to correlations between / and / + 1 multipoles, exactly as in 
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also predicts similar correlations in CMBR polarization. We analyze the CMBR temperature and polariza¬ 
tion data in order to extract the signal of these correlation between I and / + 1 multipoles. Our results for the 
case of temperature using the latest PLANCK data agree with those obtained by an earlier analysis. A de¬ 
tailed study of the correlation in the polarization data is not possible at present. Hence we restrict ourselves 
to a preliminary investigation in this case. 
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I. INTRODUCTION 


The Cosmic Microwave Background Radiation (CMB) shows a hemispherical power asymme¬ 
try, due to which the power in the two hemispheres is significantly different []Tj-8]. A dipole mod¬ 
ulation of a statistically isotropic signal provides a useful parametrization of the observed power 
asymmetry. According to the model, the observed temperature fluctuation AT along a direction n 
is expressed as [1914121. 

A? (n) = AT {n) ^1 +A%i ■ nj , (1.1) 


where AT ( h ) is a statistically isotropic field, A the dipole amplitude and X\ the dipole direction. 
Throughout this paper we shall denote the observed fields, which are assumed to have some con¬ 
tribution due to dipole modulation, with a tilde and the corresponding fields in an isotropic model 


without a tilde. Choosing our axes such that X\ is along z, Eq. (1.1) can be written as 


AT (n) — AT (w)(l+Acos0). 


( 1 . 2 ) 


As shown in [HTH T31. the two point correlation of such a modulated temperature field would show 
correlations between l and l + 1. 

If the observed signal of hemispherical anisotropy or equivalently dipole modulation is related 
to a physical effect, we expect a similar signal to be present in polarization fields too. Several 
studies have associated this effect with a primordial inhomogeneous or anisotropic model. Such 
models lead to a modification of the primordial power spectrum which culminates in depicting 
correlations between different multipoles similar to those predicted by Eq. [TT| In a recent paper 
it has been shown that such a primordial model also leads to correlations between / and / + 1 
of the polarization fields. In this paper we propose a dipole modulation model for the CMBR 
polarization field, analogous to Eq. 0 Such a model is useful to empirically characterize the 
observed hemispherical anisotropy that might be present in the polarization data, irrespective of 
the physical cause of its origin. We show that this model leads to correlations between the / 
and l + 1 multipoles for the polarization fields. We also determine the explicit form of these 
correlations. 

We search for such correlations in the recently released Planck experiment data in both tem¬ 
perature and polarization signals and compare them with previous results and predictions. In the 
case of polarization a detailed study is not possible due to difficulty in handling and interpretation 
of the noise files. Hence in this case we restrict ourselves to a preliminary investigation. 
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II. TEST FOR DIPOLE MODULATION 


The modulated temperature field is given by Eq. (1.2). This being a field on a sphere, can be 
expanded in spherical harmonics as 


A f(n) = £af m r,,„(/ 


l.m 


The two point correlation of the temperature field in multipole space can be written as 


(II. 1) 


(4,,4i> = j (*) 1/w («') (AT (n) AT (. n ')) 

As shown in lfl3l . using Eq. (1.2) we obtain 

(dj m dj,^) = CjS u ,8 mm , +A (cjt +Cf ) 

where 


eO _ c 
?lm\l'm! u mm' 


(l + m+ 1)(/ — m + 1) / + — 

(2/ + 3)(2/ + l) W + V (2/+l)(2/-l) d/V_1 


(II-2) 


(H.3) 


(H.4) 


In Eq. (II.3), the first term on RHS corresponds to the isotropic part of the correlation (d/ m d/, ; *,)i so 
and the second term is the contribution of the modulation, mo d. As we shall see the 

multipole power Cj does not get any contribution from the modulation term. This is found to be 
true also for the power in the polarization fields, to be discussed later. Hence we do not denote it 
with a tilde. We follow lfT3l and seek correlations between / and / + 1 multipoles, which can be 
expressed as, 

l~T ~T* \ A \ r T , r T -1 / (/ + m + !)(/ — m+ 1) 

( a lm a l+lm) [Q+l+Qjy (2/ + 3) (2/ + 1) ' ( ‘ 5) 

Theoretical models used to explain the dipole modulation of the temperature field predict a 
similar correlation between / and l+l multipoles in the CMB E-mode polarization field lfl4l iT5fl 
in the same direction. These predictions may be tested in future by determining these correlations 
in the CMB polarization field. 

A detailed discussion of CMB polarization is contained in [fT6lfT71l and here we use the notation 
of llT6ll . The CMB polarization field is characterized by two Stokes parameters Q and U, while the 
temperature fluctuation field corresponds to Stokes’ parameter /. Here Q and U denote the dipole 
modulated polarization fields. Under a rotation by an angle y/, the temperature field transforms as 
a scalar, while combinations of Q and U behave as spin ±2 fields on a sphere, viz. 

{Q±iU)\ii) =e* 2i V(Q±iU)(h), (II.6) 
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and can be expanded in spin ±2 harmonics as 


(Q±iU)(n) =^d± 2 ,/m ±2 Yi„ 

Im 


in). 


(II-7) 


Using the spin raising and lowering operators 5 and 9, spin 0 objects can be constructed from Q 


and U fields Ifl8l . Using 9 and 9 suitably on Eq. (II.7) we get 


9 2 (G + /f/)(n)=£ 

Im 

m-rU)(n)=Y. 

Im 


'(/ + 2)!_ y 
(/ 2)\ a2,lm Y lm' n ) 


\l + 2)L y _ 

2 y a —2,l?n lm\ n ) ■ 

Finally the standard E and B mode polarization field can be expressed as, 


(II-8) 

(II-9) 


E{n)=- 1 - [9 2 (<2 + iU) + 9 2 (Q- iU )] = £ 

~ Im 


\l + 2)\, E 


{1-2.)'. 


a,m Y lm(n) 


B(n) =- [8 2 (S+ iU) - S 2 (Q- /(/)] = £ 

Im 




( 11 . 10 ) 

(11-11) 


The coefficients af in and af m are defined as linear combinations of a± 2 jm as: 


l lm 


®lm ^ {p-2,1m T S.—2,lm) 


a B 


Im ~2 (^2,/m 3.-2,Im) 


(H-12) 

(11.13) 


The af m s define the E-mode polarization in multipole space and are unchanged under parity trans¬ 
formation in contrast to a B m s which do change sign under such a transformation. The scalar fields 
defined in Eqs. (II. 10) and ( |II.l 1| ) are the real space constructs of the E-mode and B-mode polar¬ 
izations representing the irrotational and curl components of the CMB polarizations respectively. 
In this work we are interested in the E-mode field. We define the auto correlation of the E field 
and cross correlation of the E and T fields as 


C; 


EE 




C J E =^fV\Y,( a fm S Jm) 


(11.14) 

(11.15) 


In order to study the / and l + 1 correlations we construct 


r xx _ 
K /+1 — 


1(1+1) 


m=+I 


V+VJ?-! 


E («£.«/+ iJ, 


(11.16) 
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and define our statistics as 


nXX _ 
— 


Y r xx 
L w,/+1- 


(H.17) 


Here X can be either T or E giving us TT, EE , TE and ET correlations. We search the direction 
for which the statistic S X H X maximizes in each of the maps. We also define a quantity R as the ratio 
of the anisotropic part to isotropic part, i.e., 


R 


I 

L ^max 
7 . 

1 min 


C 


xx 

/,/+! 


/(/ + i)q 


■xx 


(11.18) 


This may be seen as a measure of the fraction of the anisotropic effect to the isotropic power. 


III. DIPOLE MODULATION IN POLARIZATION 

The dipole modulated polarization fields are denoted by Q (n) and U (h) where n = (6,(j>). We 
also define 

a± («) =Q{fi)± iU (n ), (HU) 

OC±(n) — Q(n) ± iU(n) where Q and U are the standard unmodulated fields in an isotropic model, 
a — al and a ~ a* . The preferred direction X is taken to be the same for both Q and U as 
well as the temperature field lfT51 . In analogy with temperature, we propose the following model 
for dipole modulation of polarization: 

a + (n) — a + (/I) ^1 +ApX ■ 77 j , 

( 77 ) - ( 77 ) (^1 +A* P X ■ nj . (III.2) 

Here A P — A \ + 7 A 2 is a complex parameter. We choose our coordinates such that X — z and hence 
X ■ 77 = cos 0. In terms of the Stokes’ parameters, we obtain 

Q — Q(l +Ai cos 0) — UA 2 COS 0 

U - QA 2 cos 0 + 1/(1 +Ai cos 0) (III.3) 

Using Eqs. |II.7| and [II.13| for the modulated polarization fields, we obtain 

«± = -£ {4n ± i4n) ±2 Yim , (HI-4) 

l 111 

where d P [ m and denote the harmonic coefficients of the modulated fields. Inverting the above 

equation we obtain 

- {“L ± ld L) = J «± («) ±2 Yi* m (n)d£l. (m. 5 ) 
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This leads to 


af m = ~\j i a + (») (1 +A P cos 0) 2 Y,* n (n) + OC-(n)( 1 +A* P cosd)- 2 Y l * n (n)]dn : 


(III.6) 


where we have used Eq. (III.2). We also obtain a similar equation for a B 


lm ■ 


A. Correlations of the Dipole Modulated Polarization Field 


The two point correlations of the dipole modulated E field harmonic coefficients can be ex¬ 
pressed as, 

Ka £'} = \(h + h+h +h) (ni.7) 

where 

h = JJ dEldEl' (a+(n)oc-(n)') (1 +Apcos 0 ) (l -l-A^cos 0 1 ) 2 Yj* m (fi) 2 T//,„/(/ 7 / ), 
h = Jj ' d Eld El' (a + (n)a + (n )) (1 -\-Apc os 0) (l -\-Apc os 0 1 ) 2 y / * ; (n) - 2 E/' m '(« / ); 

/ 3 = JJ dCldQ! {a-(h)a-{ti)) (1 +A^cos 0 ) (l -Moose') - 2 Y* m (n) 2 T /W (n'), 

Z 4 = JJ dEldEl' (a_(n)a + (n')) (1 +Apcos 0 ) (l -l-ApcosO') 7 ) _ 2 y//„/(n / ). 

The two point correlations appearing on the right hand side of these equations can written as: 


(o+(/i)a_(*')> = E (C^ + C^) 2 y /V («) 2y; W '(«) ; 

(a + (n)a + (n')> = E (C^-C^)(-ir" 2 y /v («) 2i> ( _ w , 0 (n'), 

l"m" 

o-(*)o-(«')>= £ (c?-0(-i)”" 2 y; w ,(A') 2 y;. ( _,„»)(«), 

l"m" 

(a.(n)a + n')) = E (C^ + C^) 2 y;, (w) (n) 2^ ( _ w ,)(n')- 


l"m" 


where we have used, (af m <af* t ) = C E S u , 8 mm ,, {af m a%* Vrt/ J = Cf 8 iv 8 mrn! , (af m af m ) = 0 and 
^y^ = ( —l) m 2 y/ f _ m) . Here Cf and Cf represent the isotropic power spectrum corresponding to 
E ox B modes respectively. As we shall see the anisotropic model does not contribute to the power 
spectrum. Hence these also represent the power of the tilde fields. 


Substituting the resulting expressions of /, in Eq. III. 7 we obtain 

bn®bn 1 ) = Cf E 8 ui 8 mm r + - (lM\ +M 2 +M 3 + M 4 ) 


(III.8) 
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where M, represent the corrections due to dipole modulation and are given by, 

Mi= £ (Cff + C™) [[dado,' (Apcos 0 + A* cos 0') 2 Y lv (n) 2^ W '(«) 2 ^ w («) 2 W(«'j 

l"m" 

M 2 = £ (C^-Cf)(-1)' M " W fl d£ld&' (Apcose' + Apcosd) 2 Y Fm rin) 2 Yv'(_ m " ) (n) 

l"m" 

m 3 = £ (cf,F (,i;,c<,se+/ i;,c.k«') ,„»,(«) 

l"m" 

x 2 ^/(-,m)(”) 2 Yi' m i(h'), 

M 4 = £ (C^ + ejSSB) (-l) w ' +w H JfMn' (A^cos 0 +A/>cos 0') 2 ^ ( _ w »)(«) 2 ^ ( _ w ,)(n ) 


l"m" 




where we have assumed that the modulation parameters Ai and A? are small and dropped higher 
order terms. We can evaluate these integrals by using 

/■2;r 


./0 ./0 


2 r/m(n) 2 Y$ n t(n)dn = 8 lv 8 nm t. 


Furthermore we define 


/■27T /-;r 

I (/, m , /', m') = / / 2 F /W (n) 2 F / l m , (n) cos 0r/f2 = 5 W)W /K (/, /', m ), 

J 0 J 0 

This integral can be expressed in terms of the Wigner 3-j symbols by using 


(HI-9) 


(III. 10) 


y io(«) = cos0. 


We obtain 


(III. 11) 


I(l,m,l',m') = (-1)'"' yj (21 + 1) ( 21 ' + 1) 


/' / 1 
2-2 0 


/' / 1 

— m' m 0 


(III. 12) 


The Wigner 3-j symbol obeys the condition, 


h h h 


m\ ni 2 m 3 


= 0 if |/i-/ 2 |>/ 3 


Using this we find that 


(III. 13) 


K (/,/', m) = 0 if/' >/+l and l' < 1-1 


(III. 14) 
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For the remaining cases, l — l' and /' = / ± 1, it can be expressed as, 


K(/,/',m) = (-l) Z+Z '//(/,/',m)U (/,/', m), 


(III. 15) 


where 


) = -2 


' (2/ +1) ( 2V + 1) (/ - m )! (/ + m )! (/' - m )! (/' + m) 

(/ + 2 )! (/ — 2 )! (/' + 2 )! Jv — 2 )! 


(2Z)!(Z / + 2)!(/ / -2)!g /+ i / 2m(2/)!(/ + 2)!(/-2)!^ (2/ / )!(/-2)!(/ + 2)!g / , 1/ 

^ ’ ,m) {21 + 3)!(/ — m) !(/ + m) ! + /(/+ m)!(/-m)!(2/+ 2)! + (/ / + m)!(/ / -m)!(2/+ 1)!' 

The function IK (/,/',m) is explicitly evaluated in the next subsection. We can now write the inte¬ 
grals Mj as 


Mi — (—l) z+// 8 , 

\i+r 


,'H (/', l,m) U (/', /, m) [A P (cf, E + Cf) + A£ (cf £ + Cf B )] , 


M 2 = (-1) Z+Z SnrfAp [H (/,/',-m)U (/,/',-m) (cf £ -Cf s ) + #(/', /,m) U (/',/, m) (C^-C/F)] , 
M 3 = ( —l) z+// 8 mm rA*p [H (/,/',-m)U m) (cf E -Cf B ) m) U (/',/, m) (Cf £ -Cf s )] , 


M 4 = (-l) w 5 mra ,// (/,/', -m) U (/,/', -m) [A|> (c/F +Cf) +A/> (cf £ +Cf 5 )] . 


) z+z \ 

where we have used 


and 


U (/,/', m) = U (/',/, m) 

// (/, l',tn) = H (l ', /,m) = // (/. —m) = H (l. —m). 


This finally leads to 


1 


(W*) =Cf E 8 lv 8, nm > + 8 mm ,- [K(/,/» (ApCf £ +A;cf £ ) + K(/,/',-ro) (A^C/F+A P Cf £ )]. 

The first term on the right hand side of this equation is the standard contribution due to an 
isotropic field. The second term arises due to dipole modulation. In this term only contribu¬ 
tions linear in the dipole parameters A\ and A 2 have been kept. We see that the modulation 
model, Eq. III.2| leads to correlations between multipoles l and l + l besides also leading to 
additional contributions proportional to 8 mm i8u'. However the latter contributions cancel out after 
summing over m and can be ignored. Hence after summing over in the dipole modulation term 
leads to correlations only between l and l+l. We also notice that the terms proportional to 8 i+ 1,/' 
and 5/_i p in Eq. III.15 are symmetric under the interchange m -H- —m. Hence we deduce that 


= IK (/,/',—m). Using this we obtain 


( 5 /m 5 S) = C l $11'5mm' + A, IK(/,m) (C EE +Cf E ) . 


(III. 16) 
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where we have ignored the contributions proportional to 5; // since they cancel out after summing 
over m. Similarly for the B mode polarization, we obtain 

(5/«5/w> = Cf B 8 u ,8 mm , + 8 mm ,A^(li l',m) (cf + Cf e ) . (III.17) 

The dipole modulation model, Eq. III.2[ is very useful for characterizing a signal of anisotropic 
power that might exist in the polarization data. It allows an empirical parametrization of such 
a signal. Furthermore it can be used to perform simulations which are required for a statistical 
study of the anisotropy. We explicitly demonstrate this in the present paper. We also point out 
that an alternate model in which we may directly introduce a dipole modulation in the E mode 
polarization simply does not work. 


The correlations of the E and B mode multipoles, Eqs. III. 16 and III. 17 depend only on the 
parameter A i and are independent of A^. Hence we can directly extract A 1 by studying the E mode 
correlations. 

A similar calculation for the TE mode correlations leads to the following result 

{dj*i a E m ,) = Cj E 8n>8 mm > +AC EE tf m .j, m , + -ApCj E 8 mm iK(l,l',m) + -A* P Cf E 8 mm iK(l,l', —m ) 

(III. 18) 

In this case also, after summing over m, the S/p term in IK(/,/',m) and K(/. /'. —in) drops out. 
Hence it can be ignored and the remaining terms are symmetric under m —m. Therefore we can 
express this result as 

{^Jm^Vml) = C I E &ir&mm' + AC V tfm-l'm' + A \Cj E 8 mm iK(l,l' ,m) (III.19) 

B. Calculation of the Polarization Correlations 


We next explicitly evaluate the integral in Eq. (III. 10). The spin 2 harmonics lfl9l can be 
expressed as 


2X1,171 


, /(2;+l)(/ — »>)!(/+»>)! 

= (-1) e V 4 jt(/ + 2)!(/ —2)! 


X 




2l~2r—2+m 


cos f) 


2r+2—m 


. (III.20) 


The (j) integration in Eq. (III. 10) leads to the factor 2 7tS mm r. The 0 integral is evaluated by using 
the identity 


/o 


^CQ.asinesiWgWm -2 gp . (01.21) 


p ^ m+;i+6 ^ 
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which can be derived by using |[20ll 


r nl2 W m+l \W»+l\ 

/ sin"' 0 cos" Odd — — 2 ^ 2 

I o 2r(^±f±^) 


The function I (/, m, / , in ) defined in Eq. III. 10 can now be expressed as 


I = (—1) /+/ S mm iH (l,l',m) U (l,l',m) , 


where 


n (i ,/ mV _. (/-2)!(/ + 2)!(r + 2)!(r-2)! 

U(/,/ ,m) — 2 ( 1 '+ 1 + 2 )! S ( /,Z ’ )’ 

1 - 21'-2 

S(l,l',m) = Y, £ 

r =0 t =0 


(III. 22) 

(III.23) 
(III. 24) 


and 


& (/, l\ r. t. m) — 


( —1 ) r+t (l + l' — r — t + m — 2 )\(r + t — m + 2)\(2r + 2t — 2m + 4 — l — l') 
r\t\(l — 2 — r)\(l' — 2 — t)\(r + 2 — m)\(t + 2 — m)\(l — r + m)\(l' — t + m)\ 


We next show that for /' = /, 

§ (, l,l',m ) 

Proof: For l — /', 


§(/,/» = £ 



(III.25) 

-4m (2/)! 

m)! (/ — m)! (/ + 2)! (/ — 2)! ’ 

(m.26) 

f (-D'(,P + 2 J(T) 

' (Z-2)!(/ + 2)!(/-2)! 

(DL27) 


where 


/~ 2 /;_ 9\ I—2—r 

P=£(^1)M j (2r + 2t — 2m + 4 — 2/) x (/W+m + i)[](t-/n + 2 + v). (111.28) 

/=o V f J 


s= 1 v= 1 


We next express the two products as, 


l—2—r 


]Q (Z —f + /n + s)P[(f — m + 2 + v) = (—1) 

5=1 V =1 


1 2 r t l 2 + a\t l 3 + a2t I 4 +.. . + ai-2 


where a,- e Z. Eq. III.28 now becomes, 


1—2 


1-2 


P= («i) / - 2 -'-£(_i)M - “ ) (2r + 2f — 2m + 4 — 2/) 

r=0 ' { 


1 _9 /_T 7_4 

f “h + . . . + Cll -2 


By using Eq. VII.3 we find that only two terms, i.e. those proportional to t l 1 and t l 2 , contribute. 
Thus we obtain 

l ~ 2 t (l- 2 N 


P =( -i)'- 2 -'-£(-i)' 
t=o 


(2 r - 2m + 4 - 21 + 2a l )t , ~ 2 + 2 t l ~ l 


(III.29) 
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The constant ci\ can be determined by using the result that if 

n 

J~[(x + a,-) = V 7 + x n ~ l ai + 

i= 1 

then ci\ — I" =1 a,. Thus we obtain 


l—2—r r ^ 

ai = ( — l ~ m — s) + (—m + 2 + v) = - [—3/ 2 + 7/ — 2m/ -f 4m — 2 + 2r(2Z + 1)] . 

s=l v=l 2 

Using Eq. |VII.3| we can express Eq. |III.29| as, 


P = (-!)-'■(/ -2)! [-2/ 2 + 2/-2m/ + 2m + 4r(/ + l)] . 


Substituting in Eq. III.27 we obtain 


1—2 


§(/,/» = £ 


1 + 2 \ f l — 2\ [—2/ 2 + 21 — 2ml + 2m + 4r (/ + 1 


n\l*-r + m 
r= 0 v 1 


This sum can be divided into two parts 


-2l 2 + 2l-2ml + 2m)' x l ( 1 + 2 \ (1-2 


(l + 2 )\(l-2)\ 


(Z + 2)!(/-2)! 


I 

r=0 


l — r + m 


+ 


4(1+1) 


1-2 


1 + 2 

(l + 2) \ (l — 2)\ \l-r + m 


1 — 2 
r 


In the second sum r — 0 does not contribute. Hence after some simplifications, it can be re¬ 
expressed as, 


4 (/ — 2) (/ + 1) y? / 1 + 2 

(l + 2)1(1 — 2) \ \l — t — 1+ m 


1-3 

t 


We can evaluate both of these sums by using the Vandermonde Convolution property of binomial 
coefficients ll2Til which can be stated as, 

Tn + p\ 


I 

k=0 


P 
n — k 


n J 


, m + p>n & m, 77, p> 0. 


We finally obtain 


, (-2/ 2 + 2/-2m/+2m + 4)(2/)! 

s (/, / , m) = x , 77 : ^. /7 —+ 


4(/ + 1) (/ — 2) (2/ — 1)! 


(/ + m) ! (/ — m)! (/ + 2)! (/ — 2)! (/ + 2)! (/ - 2)! (/ - m)! (/ - 1 + m)\ 

-4m (21)1 


l(l + m)\(l—m)\(l+2)\(l — 2) \ ’ 


(111.30) 


which is the desired result and leads to U(Z, I'. m) given in Eq. III. 15 for the case / = /'. 
We next show that 


2 ( 2 /)! 


SfZjZ^m) = — < fi+m)!(/-m)!(/+2)!(/-2)! 

V ’ l 2(2/')! 

(I'+m) ! (/'—m)! (/'+2)!(/'—2)! 


/' = / + 1 
/' = /-! 


(111.31) 
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Proof: We first consider the case l' = 1 + 1. We can write Eq. (III.25) as 

1 + 2 \{l'- 2 \( 2 r+ 2 t- 2 m + A-l-l') 


= (- 1 ) 


r+t 


l — r + m 


(/'— 2)! (Z + 2)! 


x 


(l + l' — r — t + m — 2) ! 

(r + t — m + 2)! 

(/ — r — 2)! (/' — t + m)\ 

I- 

+ 

to 

1 

3^ 

l_ 


(III.32) 


After simplification of the terms in the two square brackets, this becomes 


(-D w (,'a„)(': 2 )(' , 7 2 ) 

(f'-2)!(/ + 2)!(/-2)! 


l-2-r 


(2r + 2t — 2m + A — l — l') — t+ m +s)Y\(t — m + 2 + v) 


j=l 
l -1 


V=1 


We can write the term in the square brackets above as atf 1 1 +a\t l 2 ... a/_i, where a,- G Z. Keep 


ing r fixed, the sum over t in Eq. III.24 yields, 

/'-2 V -2 n '— 0 

Y & (/,/', r,t,m) = J^(-l) ? ' 

f=0 f=0 


ClQt1 ^ + Cl\t^ ^ + ... + Gj_ i 


(III.33) 


By using Eq. |III.32| on the left hand side and by comparing both sides, we find that ciq is equal to 

we obtain 


2(— 1 ) 1 2 '. Using the second case of Eq. 


VII. 3 


l '—2 


= —2 (—I) -r (/ — 1)! = -2 (-1) r (/- 1)! 

t =0 


Finally the sum over r, after simplification, yields 


S () 


(/ + 2 )!(/ — 2 )! 


1—2 

E 

r =0 


1 + 2 
l —r + m 


1-2 


- 2 ( 2 /)! 


(/ + 2)!(/ — 2)!(/ + m)\(l — m) \ ’ 


where we have again used the Vandermonde Convolution property. A similar analysis for the case 
/' = / — 1 yields 

S (,) : 


- 2 ( 2 /') 


(/' + 2)!(/'- 2) !(/' + m)!(/'-/«)!' 


These lead to the result for U(/, l 1 . m) given in Eq. III.15 for the cases l' = l ± 1 


IV. DATA ANALYSIS 

For the case of temperature data, we perform a detailed analysis of the signal. This allows us to 
obtain updated results for the statistic Sj/ with the 2015 Planck data. We studied this statistic in 
the multipole ranges 2 < l < 64, 30 < l < 64 and 30 < l < 100 in the Planck 2015 CMB intensity 
maps. The dipole modulation signal in the CMB temperature map was observed in the multipole 
range 2 — 64 lf22ll . 
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FIG. 1. Left: SMICA inpainted temperature map. Right: SMICA temperature map with the masked portions 
tilled by isotropic randomly generated data. Both maps are for N .side = 32 and the temperature is given in 
units of K. 

For the case of CMB polarization a detailed analysis is not possible at this stage since the data 
for low l is unreliable. Furthermore we are unable to properly simulate the noise corresponding to 
PLANCK detectors. Hence, even for large /, it is not possible for us to obtain a reliable estimate of 
the errors and the significance of the signal. For this reason we confine ourselves to a preliminary 
analysis of the polarization signal. However it may still serve a useful purpose in revealing the 
preferred direction indicated by data. We searched for modulation signal in the multipole ranges 
40 — 100, 40 — 125, 50 — 100, 50 — 125, 50 — 150 and 50 — 200. The lower limit of l = 40 was 
chosen since the data for lower l is so far poorly understood. 


A. Planck 2015 temperature data analysis 

We have performed our analysis on both Commander and SMICA IQU maps. For SMICA 
maps, we have performed the analysis after masking the maps and then inpainting the masked 
maps using the MRS package of iSAP software or alternatively masking and then filling the 
masked portion of the map with isotropic data generated using the CAMB simulation package 
with Planck 2015 parameter set. The data analysis was performed with HEALPix software ll23ll . 

The analysis of the Commander and SMICA inpainted maps was performed identically. The 
maps were downgraded, without any smoothing to either A^de = 32, for multipole ranges 2 — 
64 and 30 — 64, or N s id e = 64, for 30— 100 multipole range and the analysis was performed on 
these downgraded maps. We removed the dipole and monopole from CMB intensity maps and 
calculated the quantities Y.iCff +1 and E/^(7 + 1 )Cj T over the aforementioned multipole ranges. 
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Map 

Sj/ in 10~ 2 mK 2 

A 

(l,b) 

P-value 

Commander 

2.55 ±0.68 

0.082±0.018 

(232° ± 18°, — 14° ± 18°) 

0.20% 

SMICAfinp.) 

2.39 ±0.70 

0.069 ±0.013 

(236° ±27°,-11° ±20°) 

0.70% 

SMICA(filled) 

2.44 ±0.71 

0.078 ±0.019 

(242° ±16°,-17° ±20°) 

0.50% 


TABLE I. The maximum TT Mode Sh values along with the dipole modulation parameter A, the preferred 
direction of maximization and the P-value. 


We also used another method to analyze the SMICA maps. In this case the SMICAI maps were 
first masked with their respective masks. The pixels which were masked were filled with data from 
I simulated maps generated using lensed scalar C; values generated with the CAMB Boltzmann 
solver ll24ll for the Planck 2015 parameters ll25l . and the Synfast program from the HEALPix 
package. The masked and filled I maps have A/ S jd e = 2048, same as that of the original maps. 
These maps were smoothed with a FWHM equal to 3 times the pixel size of the low resolution 
map with A s id e — 256 before downgrading the map to remove the discontinuities at the mask 
boundary. The smoothed maps are then degraded to A s j t | e = 256 followed by a further degradation 
to A^ s ide = 32 or A s ide = 64, depending on the multipole range, for final analysis. We generated 100 
such filled maps and the individual results were found to depend on the random realization used to 
fill the masked portions of the sky. The results presented here are the mean values for the statistic 
Sj/ and direction for which it maximizes. The SMICA I maps masked and reconstructed using 
these two procedures, i.e. inpainting and random filling, are shown in Fig. [T] 

We have fitted the TT mode values of Sh in order to extract the value of the dipole modulation 


amplitude A of Eq. (1.2). We simulated 100 isotropic CMB maps using Planck 2015 parameters 
with A s ide = 512 (A s ide = 1024 for SMICA filled analysis). These maps were rotated to have the 
z axis pointing along the direction of maximum statistic and they were modulated using Eq. CD- 
The modulated maps were downgraded to A s id e = 32. Each of the downgraded simulated maps 
were fitted for the value of A that would give the value of Sh closest to the one observed in the data 
along the direction along which it maximizes in the Planck 2015 CMB maps. We averaged over 
100 best fit values of A obtained by this method giving the modulation amplitude for the results 
given in Table [IJ 

To estimate the error in S ! H r we generated 1000 maps at A s id e = 512 (A s id e = 1024 for SMICA 
filled analysis) and modulated the maps with the best-fit value of A along the direction of maxi- 
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I Range 

S t h t in 10 2 mK 2 

A 

V,b) 

P-value 

R 

2-64 

2.55 ±0.68 

0.082±0.018 

(232° ±18°,-14° ±18°) 

0.20% 

0.065 

30-64 

1.00 ±0.43 

0.052 ±0.019 

(194° ±25°,-4° ±24°) 

66.6% 

0.040 

30-100 

0.91 ±0.72 

0.018±0.011 

(277° ±81°,4° ±30°) 

83.0% 

0.013 


TABLE II. TT mode Sh results for Planck Commander 2015 maps in different multipole ranges. 


mum statistic using relation ( |I.2| ). The modulated CMB maps were downgraded to ACde = 32 or 
A^side = 64 depending on the multipole range under consideration and Sj/ was calculated along 
the direction of modulation. The standard deviation of the 1000 simulated maps gives the error 
in Sj/. For the masked and filled SMICA maps, the process of filling the masked region of the 
SMICA maps with isotopic data is expected to introduce bias in the obtained results for Sh- This 
bias correction is relatively small l(22ll and we ignore it in our analysis. It is expected to enhance 
the signal by about 8%. 

The error estimation in the preferred direction was performed by simulating 50 isotropic tem¬ 
perature maps. We modulated these maps with the best-fit value of A along the observed direction 
in the data. The direction along which the statistic S ! H r maximises in the maps was found and the 
standard deviation of the coordinates gave corresponding errors. 

Finally to test the significance of our results we simulate 2000 isotropic CMB maps for 2 — 64 
multipole range and 500 maps for the other multipole ranges, with Planck 2015 parameters and 
search for the direction along which the statistic maximizes. The values of Sh obtained by this 
process is used to obtain probability distribution of the statistic Sh for the isotropic hypothesis. 
The histogram is given in Fig. [2] P-values for individual results were obtained as the percentage 
of simulation results that equal or exceed the observed result for the statistic. 


B. Planck 2015 polarization data analysis 

For the case of polarization, as explained above, we are unable to perform a detailed analysis 
of dipole modulation. Here we confine ourselves to simply making an estimate of the statistic Sff 
and the corresponding preferred direction using the Commander map. Since the low / multipoles 
are not expected to be reliable we confine our study to the multipole ranges 40 — 100, 40 — 125, 
50 — 100 and 50—125. Furthermore we use the dipole modulation model for polarization in order 
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FIG. 2. Histogram of TT Sh values obtained by using isotropic ACDM simulations. The observed Sh 
values are indicated by dashed lines. The results obtained by using the Commander map (blue line), SMICA 
inpainted (red line) and SMICA filled map (green line) are shown. 


to generate simulated maps which display polarization power anisotropy and to determine the 
distributions of the corresponding statistic Sh for the E mode polarization. The main purpose of 
this study is to illustrate the utility of this model. 


V. RESULTS 

In this section we first present the results for the temperature analysis and later those of polar¬ 
ization. 


A. Temperature 

The TT mode results for multipole range 2 < / < 64 are summarized in Table [I| As stated 
in the previous section the dipole modulation in the CMB temperature signal is present in lower 
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FIG. 3. From the left Sj/ sky maps for Commander data in multipole ranges 2 — 64, 30 — 64 and 30 — 100. 

multipoles up to l = 64. Our primary TT mode results are for this multipole range. We can 
compare our results with [[22| to look for changes between the Planck 2013 and Planck 2015 data. 
The results for SMICA filled maps for Planck 2013 data are: Sj { T — (2.1 ±0.5) x 10~ 2 mK 2 , a bias 
corrected value of (2.3±0.6) x 10~ 2 mK 2 , along (229°, — 16°) with A = 0.074±0.019. Comparing 
the 2013 and 2015 SMICA filled map results we notice agreement in both the values of Sj/ and 
direction of maximization, while the value of A is also comparable within the error limits. The 
results for SMICA inpainted map for Planck 2013 data are: S ! H r — (2.7 ±0.7) x 10 2 mK 2 . along 
(232°,—12°). For SMICA inpainted maps too, we find good agreement with previous results. 
We also note that the 2015 results have smaller P-values when compared with isotropic ACDM 
simulations generated using 2015 Planck parameters. We can compare our results with Planck 
Collaboration’s analysis of dipole modulation [|26l . Planck 2015 best fit values of modulation 
amplitude A is 0.063 ^[|'q 2 3 for Commander maps and 0.062 ^qq 2 3 for SMICA map. The direction 
of modulation is (213°, —26°) ± 28° for both SMICA and Commander maps. Both the modulation 
amplitude and the direction agree with our results within the quoted errors. Our results are also 
consistent with those obtained in ll27l . 

It has been found by an earlier analysis that the hemispherical anisotropy effect rapidly dies 
out when the lower multipoles are excluded or when summed to higher multipoles [261 1281 - 43011 . 
In order to study the multipole dependence we also investigate the effect in the multipole ranges 
30 — 64 and 30 — 100. The S\ { T obtained in 30-64 for Commander map is contained in Table [nj 
In this case we do not find a significant signal of 1,1 + 1 correlation. The temperature map Sh 
is very much in agreement with expectations from isotropic theory. We also notice that the S l H r 
maximizes along directions which change from one range to another. This can also be seen from 
Fig. [ 3 ] These maps show S ! H r for different directions in ranges 2-64, 30-64 and 30-100. It can be 
seen that left to right the dipole pattern becomes less distinct. 

The SMICA inpainted TT mode results are as follows: 
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• In the range 30-64, Sj/ maximizes at (1.24 ±0.43) x 10 2 mK 2 along (191° ± 27°, —4° ± 
28°) with a P-value of 45% and A = 0.059 ± 0.019. 

• For 30-100 range, Sj/ maximizes at (1.45 ±0.71) x 10 2 mK 2 along (207° ±38°, —18° ± 
27°) with a P-value of 41.2% and A = 0.027 ± 0.013. 

The SMICA filled map temperature results are: 

• In 30-64 range, Sj/ maximizes at (1.00 ±0.43) x 10 2 mK 2 along (228° ± 28°, —9° ±27°) 
with a P-value of 66% and A = 0.055 ±0.021. 

• The 30-100 range Sj/ maximizes at (1.64 ±0.77) x 10~ 2 mK 2 along (238° ± 34°,—11° ± 
28°) with a P-value of 27% and A = 0.031 ±0.014. 

Both SMICA inpainted and filled maps results show similar patterns as Commander map results. 


B. Polarization 


In this section we present our results for polarization. As explained above, we do not make an 
attempt to compute either the errors or the significance of this effect due to uncertainties in the 
noise simulation in polarization. However we do present the results of a simulation in order to 


illustrate the utility of the polarization dipole modulation model, Eq. III.2 

we give the results for Sff, the preferred direction and the ratio R for the E mode 


In Table 


III 


polarization. The result for the case of the multipole range 40 — 100 is mildly interesting. This 
is because the preferred direction aligns closely with the CMBR dipole. A preferred direction 
similar to the one obtained in the range 40 — 100 has been seen in many other studies OTl 32], 
including the radio polarization dipole axis 0511 . the CMB quadrupole - octopole alignment axis 
031 , the NVSS dipole [|35T- l38l . the radio polarization flux dipole axis ||39| as well as the optical 
polarizations from distant quasars |3~1.140l - l42l . It is possible that the present signal in E mode 
polarization in this range arises from some residual contamination of the low l noise systematic 
bias. Alternatively it may be a signal of some astrophysical or cosmological effect. This may be 
settled by future refinements in data. 

For the higher multipole ranges the preferred direction starts to deviate. For example, in the 
multipole range 40 — 125, 50 — 100, 50— 125 it lies closer to the galactic plane. This might be 
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I Range 

S% E in 10 6 mK 2 

(l,b) 

R 

40-100 

6.6 

(260°,44°) 

0.031 

40-125 

9.4 

(291°,14°) 

0.023 

50-100 

6.6 

(286°, 15°) 

0.033 

50-125 

9.5 

(291°, 14°) 

0.025 


TABLE III. EE mode Sh results for Planck Commander 2015 maps in different multipole ranges. 



FIG. 4. The Sff sky maps for Commander data in multipole ranges 40 — 100 (left) and 40 — 125 (right). 

an indication that it is moving closer to the axis obtained in the case of temperature. Alterna¬ 
tively, since it lies very close to the galactic plane, the signal in this range might be dominated by 
foregrounds. The plots of Sff for the multipole range 40—100 and 40—125 are shown in Fig. [4j 
We next generate 1000 simulated Q and U mode polarization maps which include dipole mod¬ 
ulation with parameter A\ = 0.05 for the multipole range 40 — 100. We compute the statistics S l fj‘ 
for these simulated maps. The resulting distribution of this statistics is shown in Fig. [5j We find 
that the distribution is close to normal. The parameter chosen is not too far from what might be 
required to fit the E mode results given in Table [Till 


VI. CONCLUSIONS 

This work has updated the previous results for dipole modulation in CMB temperature field. 
We find that the dipole modulation of the temperature field, observed in WMAP data and Planck 
2013 data is also present in Planck 2015 data. We have computed the quantity Sj/ and the results 
are found to be in agreement with the values obtained from the previous data lf22ll . Our result for 
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FIG. 5. The Sff histogram for the simulated dipole modulated E mode polarization maps with the parameter 
A\ =0.05. 

modulation parameter A and modulation direction are compatible with Planck 2015 results lf26ll . 

We have presented a preliminary analysis of the dipole modulation in the CMBR polarization. 
In this case we coniine our analysis only to the multipole values / > 40 since the lower / values are 
not expected to be reliable. The detector noise is expected to make considerable contribution in this 
case. Since this information is not available to us, we do not attempt to compute the significance of 
the signal or the associated error in extracted parameters. In any case our preliminary investigation 
reveals some interesting results about the direction of the signal. We find that for the range of 
multipoles 40—100 the preferred direction is close to the CMBR dipole and for a slightly higher 
range 50 — 100 or 40 — 125, the direction shifts closer to the galactic plane. The direction found 
for the lower range of multipoles is somewhat interesting since a similar direction has been found 
in several other data sets OTI -[39ll. 

For the case of polarization we also propose a dipole modulation model, Eq. |III.2| This is a 
useful model which generalizes the temperature dipole modulation [9l4l2|. Eq. 0 to the case 
of polarization. We find that this model leads to correlations between / and l + 1 multipoles of 
the polarization field. We determine the form of these correlations for the case of EE, BB and 
TE fields. We show the utility of this model by creating simulated polarization maps with non- 
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zero values of the dipole modulation parameter, A\. We find that distribution of the resulting 
statistic for the case of E model polarization is well described by a Gaussian. We expect that this 
dipole modulation model and our proposed analysis procedure may be very useful in the study of 
anisotropy which might exist in the CMBR polarized field. 


VII. APPENDIX 


In this Appendix we derive the mathematical results required in proving the results in section 


III B Let f(n,p) be defined as 


f(n,p) = JV’(-l)' 

r =0 

where n. p > 0. This function satisfies the recurrence relation 

f( , ,, , (-")Ej=o f(n~l,q){ p ) n^O 


(VII. 1) 


(VII. 2) 


0 n = 0 

Proof: The result for n = 0 can be verified by direct substitution. For n > 0, we can write f{n. p+\) 


as 


Y r P. 


(77-1)! 


/(„./,+ !)=„ E(-l)'. {n _ r)x(r _ l)V 

On the right hand side the term corresponding to r — 0 is zero due to the factor (r — 1)! in the 
denominator. Hence we can start the sum from r = 1. Setting r — 1 = t, we obtain 

V 

q= " 


n— 1 


/(77,p + l) = -n £(-l) f (l+ry 
7=0 


(n — 1)! A (n-\ 




t q 

1 5 


where we have used the binomial theorem. Furthermore 

f(n-l,q) = J^(-l) 


7 / n ~ 1 , tq 


7=0 


and hence 


f(n,p+l) = -n J^f(n-l,q) 
? =0 

This proves the result in Eq. |VII.2 for n / 0. 


The function f(n,p) defined in Eq. VII. 1 is given by 

( 


f(n,p) = < 


0 p < 77, 77 ^ 0 

(-1 ) ,! 77! p = n, 77 >0 

) n n\ p = n + \,n>0 


(VII.3) 
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Proof: The first part can be proven by using the Corollary 2 of [43] with x = 0 and p — n — j. The 
result follows since p < n when 1 < j < n. Furthermore the second part can be obtained by using 
Theorem 1 of [[431 with x = 0. 


Finally we consider the last part. By using Eq. VII.2 we find that for n — 0, p — n + 1, 


f(0,p ) = 0, which agrees with the result given in Eq. VII.3 We next consider n > 0. By using 


Eq. |VII.2| we obtain 

f(n,n+l) = -n^f(n-l,q) 

<7=0 

That is, the only nonvanishing terms in this sum are obtained by setting q = n — 1 and q — n. Using 


= —77 


/(n-l,n-l)(^_ 1 ) +f(n-l,n) 


the second case of Eq. VII.3 this can be further simplified in the form of the following recurrence 
relation: 


P(n) — —77 [— (—1) 77 ! +P(n — 1)] , 


(VII.4) 


where P(ri) — f(n,n + 1). We next proceed by induction. For n = 1, by using VII. 1 we obtain 


P(l)=/(1,2) = -1, 


which agrees with the result given in Eq. VII.3 We next assume that this result is true for n = k, 

P{k) = MET! (_!)»*! 


and show that it is also true for P (k + 1). By using recurrence relation VII.4 we obtain 

k(k+ 1) 


P(ik+l) = -(ik+l) 


-(-1 ) k+l {k+ 1)! + K ^ \ -l) k kl 


(fc+ 1 Kfc+2) ( _i )t+l(fc + 1)! 


which agrees with Eq. VII.3 Hence the third part is also proven by induction. 
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